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Abstract 



In this paper, the method of approximate transformation groups which was proposed by Baikov, 
Gazizov and Ibragimov [TJ [2J, is extended on Hamiltonian and bi-Hamiltonian systems of evolution 
equations. Indeed, as a main consequence, this extended procedure is applied in order to compute 
the approximate conservation laws and approximate recursion operators corresponding to these type 
of equations. In particular, as an application, a comprehensive analysis of the problem of approximate 
conservation laws and approximate recursion operators associated to the Gardner equation with the small 
parameters is presented. 
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(N ! 1 Introduction 

The investigation of the exact solutions of nonlinear evolution equations has a fundamental role in the non- 
linear physical phenomena. One of the significant and systematic methods for obtaining special solutions of 
systems of nonlinear differential equations is the classical symmetries method, also called Lie group analy- 
sis. This well known approach originated at the end of nineteenth century from the pioneering work of Lie 
[7]. The fact that symmetry reductions for many PDEs can not be determined via the classical symmetry 
method, motivated the creation of several generalizations of the classical Lie group approach for symmetry 
reductions. Consequently, several alternative reduction methods have been proposed, going beyond Lie's 
classical procedure and providing further solutions. One of these techniques which is extremely applied par- 
ticularly for nonlinear problems is perturbation analysis. It is worth mentioning that sometimes differential 
equations which appear in mathematical modelings are presented with terms involving a parameter called 
the perturbed term. Because of the instability of the Lie point symmetries with respect to perturbation of 
coefficients of differential equations, a new class of symmetries has been created for such equations, which are 
known as approximate (perturbed) symmetries. In the last century, in order to have the utmost result from 
the methods, combination of Lie symmetry method and perturbations are investigated and two different 
so called approximate symmetry methods(ASM) have been developed. The first method is due to Baikov, 
Gazizov and Ibragimov [TJ [5] . The second procedure was proposed by Fushchich and Shtelen [5] and later 
followed by Euler et al [3j [4] . This method is generally based on the perturbation of dependent variables. In 
[101 111) , a comprehensive comparison of these two methods is presented. 

As it is well known, Hamiltonian systems of differential equations are one of the famous and significant 
concepts in physics. These important systems appear in the various fields of physics such as motion of rigid 
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bodies, celestial mechanics, quantization theory, fluid mechanics, plasma physics, etc. Due to the significance 
of Hamiltonian structures, in this paper, by applying the linear behavior of the Euler operator, character- 
istics, prolongation and Frechet derivative of vector fields, we have extended ASM on the Hamiltonian and 
bi-Hamiltonian systems of evolution equations, in order to investigate the interplay between approximate 
symmetry groups, approximate conservation laws and approximate recursion operators. 
The structure of the present paper is as follows: In section 2, some necessary preliminaries regarding to 
the Hamiltonian structures are presented. In section 3, a comprehensive investigation of the approximate 
Hamiltonian symmetries and approximate conservation laws associated to the perturbed evolution equations 
is proposed. Also, as an application of this procedure, approximate Hamiltonian symmetry groups, approxi- 
mate bi-Hamiltonian Structures and approximate conservation laws of the Gardner equation are computed. 
In section 4, the approximate recursion operators are studied and the proposed technique is implemented 
for the Gardner equation as an application. Finally, some concluding remarks are mentioned at the end of 
the paper. 



2 Preliminaries 

In this section, we will mention some necessary preliminaries regarding to Hamiltonian structures. In order 
to be familiar with the general concepts of the ASM, refer to [BJ. It is also worth mentioning that most of 
this paper's definitions, theorems and techniques regarding to Hamiltonian and bi-Hamiltonian structures 
are inspired from [5]. 

Let M a X x U denote a fixed connected open subset of the space of independent and dependent variables 
x = (x 1 , • • • ,x p ) and u = (it 1 , • • ■ ,u q ). The algebra of differential functions P(x,u^) = p[u] over M is 
denoted by srf . We further define ,c/' to be the vector space of /-tuples of differential functions, P[u] = 
(Pi [it], . . . , Pi [it]), where each Pj e srf . 

A generalized vector field will be a (formal) expression of the form 

l=l a = l 

in which £ z and <p a are smooth differential functions. The Prolonged generalized vector field can be defined 
as follows: 

pr (n) V = V+ J J] KWl^, 



a = l ttJsSn Ua J 



whose coefficients are determined by the formula 



D J U a -j]eut) + f i eu%, (2.2) 



i=l i=l 



with the same notation as before. Given a generalized vector field v, its infinite prolongation (or briefly 
prolongation) is the formally infinite sum 

Pr"£^+£l>^, (2-3) 

2=1 0=1 J J 

where each (f) J a is given by (|2.2j) . and the sum in (|2.3[) now extends over all multi- indices J = (ji, ■ ■ ■ , jk) for 

k 55 o, l s; j k p. 

A generalized vector field v is a generalized infinitesimal symmetry of a system of differential equations 



A v [u] = A u (x,u^) = 0, u = l,...,l, 
if and only if 

prv[A„] = 0, v = 1,...,/, (2.4) 

for every smooth solution u = f{x). 
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Among all the generalized vector fields, those in which the coefficients C[ u ] °f the are zero play a 
distinguished role. Let Q[u] = (Qi[u], . . . ,Q q [u]) e stf q be a g-tuple of differential functions. The generalized 
vector field 

v Q = 2 0«M^ 

a— 1 

is called an evolutionary vector field, and Q is called its characteristic. 

A manifold M with a Poisson bracket is called a Poisson manifold, the bracket defining a Poisson structure 
on M. Let M be a Poisson manifold and H : M — > R be a smooth function. The Hamiltonian vector field 
associated with H is the unique smooth vector field vh on M satisfying the following identity 

v u = {F,H} = -{H,F} (2.5) 

for every smooth function F : M — *■ R. The equations governing the flow of vh are referred to as Hamilton's 
equations for the "Hamiltonian" function H. 

Let x = (x 1 , . . . , x m ) be local coordinates on M and H(x) be a real- valued function. The following basic 
formula can be obtained for the Poisson bracket. 

mm ATT 

In other words, in order to compute the Poisson bracket of any pair of functions in some given set of local 
coordinates, it suffices to know the Poisson brackets between the coordinate functions themselves. These 
basic brackets, 

J»(x) = {x i ,xi}, i,j = l,...,m, (2.7) 

are called the structure functions of the Poisson manifold M relative to the given local coordinates, and serve 
to uniquely determine the Poisson structure itself. For convenience, we assemble the structure functions into 
a skew- symmetric m x m matrix J(x), called the structure matrix of M. Using V-ff to denote the (column) 
gradient vector for H , the local coordinate form (I2.6[) for the Poisson bracket can be written as 

{F, H} = VF ■ J\7H. (2.8) 

Therefore, in the given coordinate chart, Hamilton's equations take the form 

^ = J{x)VH{x). (2.9) 
Alternatively, using (12.61) . we could write this in the "bracket form" as follows: 

dx i m 

the i-th component of the right-hand side being {x l , H}. A system of first order ordinary differential equations 
is said to be a Hamiltonian system if there is a Hamiltonian function H(x) and a matrix of functions J(x) 
determining a Poisson bracket (|2.8p whereby the system takes the form (|2.9j) . 
If 

j 

is a differential operator, its (formal) adjoint is the differential operator 29* which satisfies 



f P-S'Qdx= \ Q ■ @*Pdx 
Jo Jn 



)n Jn 

for every pair of differential functions P,Qe&/ which vanish when u = 0. Also, for every domain fi c R p 
and every function u = f(x) of compact support in f2. An operator 3> is self-adjoint if 2* = 2; it is 
skew-adjoint if 2>* = — 2>. 

The principal innovations needed to convert a Hamiltonian system of ordinary differential equations (|2.9[) 
to a Hamiltonian system of evolution equations are as follows (refer to 9 for more details): 
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(i) replacing the Hamiltonian function H(x) by a Hamiltonian functional Jtf?\u\, 

(ii) replacing the vector gradient operation Vif by the variational derivative 8Jf? of the Hamiltonian 
functional, and 

(iii) replacing the skew-symmetric matrix J(x) by a skew-adjoint differential operator S> which may depend 
on u. 

The resulting Hamiltonian system will take the form 

at 

Clearly, for a candidate Hamiltonian operator @ the correct expression for the corresponding Poison bracket 
has the form 

{^>,JS?} = J 8&-®8&dx, (2.10) 

whenever 3 s , S£ e & are functionals. Off course, the Hamiltonian operator ^ must satisfy certain further 
restrictions in order that (|2.10p be a true Poisson bracket. A linear operator 3 : sd q — » srf q is called 
Hamiltonian if its Poisson bracket (|2.10| satisfies the conditions of skew-Symmetry and the Jacobi identity. 

Proposition 2.1. Let <$) be a Hamiltonian operator with Poisson bracket (|2.10l) . To each functional = 
\ H dx G J£", there is an evolutionary vector field pryjf , called the Hamiltonian vector field associated with 
Jtif, which for all functionals 3> e & satisfies the following identity: 

prv#(^) = {&, JV} (2.11) 

Indeed, has characteristic '38M' = &E(H), in which E is Euler operator. (Proposition 7.2 of [5]) 



3 Approximate Hamiltonian Symmetries and Approximate 
Conservation Laws 

Consider a system of perturbed evolution equations 

^=P[u,e] (3.1) 

in which P[u, e] = P{x,vS n \e) e stf q , x e W\ ueW and e is a parameter. 

Substituting according to (|3.ip and its derivatives, we see that any evolutionary symmetry must be 
equivalent to one whose characteristic Q[u,e] = Q(x,t,u^ m \e) depends only on x, t, u, e and the x- 
derivatives of u. On the other hand, (|3.1[) itself can be considered as the equations corresponding to the flow 
exp(^Vp) of the evolutionary vector field with characteristic P. The symmetry criterion (I2.4p . which in this 
case is 

DtQv = pr Vq(P w ) + o(e p ), v = l,...,q, (3.2) 

can be readily seen to be equivalent to the following Lie bracket condition on the two approximate generalized 
vector fields. Indeed, this point generalizes the correspondence between symmetries of systems of first order 
perturbed ordinary differential equations and the Lie bracket of the corresponding vector fields. 

Considering above assumptions, some useful relevant theorems and definitions could be rewritten as 
follows: 

Proposition 3.1. An approximate evolutionary vector field vq is a symmetry of the system of perturbed 
evolution equations u t = P[u, e] if and only if 

^ + [vp,vg]=o( £ ") (3.3) 

holds identically in (x, t, u^ m \ e). (Here dvQ/dt denotes the evolutionary vector field with characteristic 

dQ/dt.) 
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Any approximate conservation law of a system of perturbed evolution equations takes the form 

D t T + DivX = o(e p ), (3.4) 

in which Div denotes spatial divergence. Without loss of generality, the conserved density T(x,t, u^ n \e) 
can be assumed to depend only on x-derivatives of u. Equivalently, for n a X , the functional 



&~[t;u,e] = T(x,t,u {n \e)dx 
Jn 



is a constant, independent of t, for all solutions u such that T(x, t, u^ n \ e) — * as x — * dfl. Note that 
if T(x,t,u( n \e) is any such differential function, and u is a solution of the perturbed evolutionary system 
u t = P\u, e], then 

D t T*d t T + prv p (T), 

where dt = d/dt denotes the partial t-derivative. Thus T is the density for a conservation law of the system 
if and only if its associated functional & satisfies the following identity 

dSTjdt + prv p (5) = o(s p ). (3.5) 

In the case that our system is of Hamiltonian form, the bracket relation (|2.1ip immediately leads to the 
Noether relation between approximate Hamiltonian symmetries and approximate conservation laws. 

Definition 3.2. Let 2$ be a q x q approximate Hamiltonian differential operator. An approximate distin- 
guished functional for 9 is a functional Sf e & satisfying £F<5Sf = o(e p ) for all x, u. 

In other words, the Hamiltonian system corresponding to a distinguished functional is completely trivial: 
u t = 0. 

Now, according to [5], the perturbed Hamiltonian version of Noether 's theorem can be presented as 
follows: 

Theorem 3.3. Let u t = be a Hamiltonian system of perturbed evolution equations. An approximate 

Hamiltonian vector field v^> with characteristic Qib~2P,2? e J^, determines an approximate generalized 
symmetry group of the system if and only if there is an equivalent functional 3? » £P — 'S differing only from 
by a time-dependent approximate distinguished functional ^[t; it, e], such that determines an approximate 
conservation law. 



Example 3.4. The Gardner equation 



u t = 6(u + eu 



2 N 



m, - u 



XXX • 



can in fact be written in Hamiltonian form in two distinct ways. Firstly, we see 

u t = D x (iu 2 + 2eu 3 - u xx ) = 9>b3#{, 

where & = D. x and 

2 



3K[u,e] = ^(u 3 + £ -u 4 + ^)dx 



is an approximate conservation law. Note that @ is certainly skew-adjoint, and Hamiltonian. The Poisson 
bracket is 

{^>,^} = f ■ D x (S^)dx 

The second Hamiltonian form is 

u t = (AuD x + 2u x + 3e(uu x + u 2 D x ) - D^ju = S8J^ Q , 



in which 



J#o [u, e] = J" — u 2 dx 
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$ is skew-adjoint and satisfies the Jacobi identity. Therefore it is Hamiltonian. 

In [8] , we have comprehensively analyzed the problem of approximate symmetries for the Gardner equa- 
tion. We have shown that the approximate symmetries of the Gardner equation are given by the following 
generators: 

vi = 8 X , v 2 = d t , v 3 = 6td x + (2eu - l)d u , 
v 4 = evi, v 5 = ev 2 , v 6 = s(6td x — d u ) = ev 3 , v 7 = £(xd x + 3td t - 2ud u ), 
with corresponding characteristics 

Qi = u x , Q 2 = 6 (it + eu 2 )u x - u xxx , Q 3 = 6tu x + 1 - 2eu 
Qi = eQi = eu x , Q 5 = sQ 2 = e(6uu x - u xxx ), Q 6 = £<?3 = e(6tu x + 1) 
Q 7 = e \2u + xu x + 3t(6uu x - u xxx )j , 

(up to sign). 

For the first Hamiltonian operator 3) = D x , there is one independent nontrivial approximate distinguished 
functional, the mass £^0 = M = \ udx which is consequently approximately conserved. 
For the above seven characteristics, we have 

QtaDzSPi, i = 1,2,4,5,6, (3.6) 

with the following approximately conserved functionals: 



u 2 dx, 



9>x = J%[u, e] = J \x 2 dx, &> 2 = -M[u, e] = J (u 3 + |u 4 + -it 2 ^ dx, ^ 4 = e^i = J | 

3*5 = e^2 = e J (v? + ^ulj dx, ^ 6 = e J (3tu 2 + xu) dx. 

For the second Hamiltonian operator S = AuD x + 2u x + 3e(uu x + u 2 D x ) — D x , 

Qi^gS&i, 1 = 2,4,5,7, (3.7) 
the following approximately conserved functionals are the corresponding approximate conservation laws: 

&2 = &i=\^\u 2 dx, # 4 = = | judx, 

0> 5 = e & 2 = e\ j \u 2 dx, & 7 = = | J (Mu 2 + xu) dx. 

In this case, nothing new is obtained. Note that the other approximate conservation law =^5 did not arise from 
one of the geometrical symmetries. According to Theorem l3.31 however, there is an approximate Hamiltonian 
symmetry which gives rises to it, namely v^> 5 . The characteristic of this approximate generalized symmetry 

is 

Q 5 rb dfS&'s = Se{iu 2 - u xx ) ss s(u xxxxx - I0uu xxx - 2Qu x u xx + 30u 2 w^J 
Note that Q5 happens to satisfy the Hamiltonian condition (|3.6|l for @ with the following functional 



u xx ~ 5 u 2 u xx + 5 u )dx 

Consequently, another approximate conservation law is provided for the Gardner equation. 

Keeping on this procedure recursively, further approximate conservation laws could be generated. But, 
this procedure will be done in the next section by applying approximate recursion operators. 
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4 Approximate Recursion Operators 

Definition 4.1. Let A be a system of perturbed differential equations. An approximate recursion operator 
for A is a linear operator 3H : srf q — > sd q in the space of g-tuples of differential functions with the property 
that whenever vq is an approximate evolutionary symmetry of A, so Vg is with Q «s S?Q. 

For nonlinear perturbed systems, there is an analogous criterion for a differential operator to be an 
approximate recursion operator, but to state it we need to introduce the notion of the (formal) Frechet 
derivative of a differential function. 

Definition 4.2. Let P[u,e] = P(x, «'"',£) e s/ r be an r-tuple of differential functions. The Frechet 
derivative of P is the perturbed differential operator D p : srf q —* stf r defined so that 

T>p(Q) = 4- P[u + eQ[u,e]] (4.1) 
ae e=o 

for any Q e g/ q . 

Proposition 4.3. If P e srf r and Q e srf q then 

D P (Q)*prv Q (P). (4.2) 

Theorem 4.4. Suppose that A[u,e] = be a system of perturbed differential equations. If M : srf q —> stf q 
is a linear operator such that for all solutions u of A, 

D A 31! ■ D A (4.3) 

where 3% : si q — > si q is a linear differential operator, then is an approximate recursion operator for the 
system. 

Suppose that A[u, s] = u t — K[u,e] is a perturbed evolution equation. Then Da = D t — Da-. If 3$ is 
an approximate recursion operator, then it is not hard to observe that the operator 3? in (|4.3p must be the 
same as Si. Therefore, the condition (|4.3p in this case reduces to the commutator condition 

St*\p K ,a\ (4.4) 

for an approximate recursion operator of a perturbed evolution equation. 

From (|4.4[) . we can conclude that if Si is an approximate recursion operator, then for all Z ^ 1 in which 
E l 3$ # 0, e l S£ is an approximate recursion operator 

(e l M) t = e l M t « e l [D K ,3t] « [T> K ,s l M] . 

In order to illustrate the significance of the above theorem, we discuss a couple of examples, including 
the potential Burgers' equation and the Gardner equation. In the first example, we apply some technical 
methods, used in Examples 5.8 and 5.30 of [9]. 

Example 4.5. Consider the potential Burgers' equation 

u t = u xx + eu\. 

As mentioned in [10] . approximate symmetries of the potential Burgers' equation are given by the following 
twelve vector fields 

2 2 
U U 

vi = d x , v 2 = d t , v 3 = xd x + 2td t , v 4 = 2td x - (xu - et— )8 U , v 5 = [u - et—)8 u , 

u 2 

v 6 = Axtd x + 4t 2 d t - (x 2 + 2t)(u - et—)d u , v 7 = evi, v 8 = ev 2 , v 9 = e(xd x + 2td t ) = ev 3 , 
vio = e(2td x — xud u ) = £V4, vn = eud u = evs, V12 = e(Axtd x + 4t 2 dt — (x 2 + 2t)ud u ) = evq, 
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plus the infinite family of vector fields 



f(x,t)(l - eu) + eg{x,t)\d u , 



where /, g are arbitrary solutions of the heat equation u t = u xx . 

The corresponding characteristics of the first twelve approximate symmetries are 

u 2 

Qi = u x , Qi = u xx + eu x , Q 3 = xu x + 2t{u xx + su x ), Q 4 = xu + 2tu x — et— , 

i i 
u u 
Q 5 =u-et—, Q 6 = (x 2 + 2t)(u - et—) + 4xtu x + 4t 2 (u xx + eul), Q 7 = eQ\ = eu x , 

Qs = eQ2 = su xx , Q 9 = eQ 3 = e(xu x + 2tu xx ), Q w = eQ 4 = e(xu + 2tu x ), 
Qn = £(?5 = eu, Q12 = sQ 6 = e((x 2 + 2t)u + 4xtu x + 4t 2 u xx ), 

(up to sign). 

Inspection of Q\, Qi, Q7, Q% leads us to the conjecture that M\ = D x + eu x is an approximate recursion 
operator, since Q3 = £%\Q\, Qg = fflxQ-j, etc. To prove this, we note that the Frechet derivative for the 
right-hand side of potential Burgers' equation is 

D K = Dl + 2eu x D x . 

We must verify (|4.4p . The time derivative of the first approximate recursion operator 3$\ on the solutions of 
the potential Burgers' equation is the multiplication operator 

(^i)t = (D x + eu x ) t = eu xt = e(u xxx + 2eu x u xx ) = eu xxx . 

On the other hand, the commutator is computed using Leibniz' rule for differential operators: 

[D K ,Mx] = eu xxx . 

Comparing these two verifies (|4.4p and proves that £%\ is an approximate recursion operator for the potential 
Burgers' equation. 

There is thus an infinite hierarchy of approximate symmetries, with characteristics 2#\Q\, k = 0, 1, 2, . . . 
For example, the next characteristic after Q12 in the sequence is 



1Q12 = e((z 2 + 6t)u x + 2x(u + 2tu xx ) + 4t 2 u xxx ^ . 



To obtain the characteristics depending on x and t, we require a second approximate recursion operator, 
which by inspection, we guess to be 

^2 = tM x + ~. 
Using the fact that ^1 satisfies (|4.4[) . we readily find 

{M 2 )t = t{Mx)t+Mx = t\D K ,M{\+® x , 

whereas 

[D K ,M 2 ] = t{D K ,M{\ + [D 2 X + 2eu x D x , ^x] 

= t[D K ,Mi] + {D x +eu x )= t[D Kl M 1 }+@ 1 , 

proving that ^2 is also an approximate recursion operator. There is thus a doubly infinite hierarchy of 
approximate generalized symmetries of potential Burgers' equation, with characteristics M l 2 ^%iQi, k,l > 0. 
For instance, Q2 = £%\Q\, Q3 = 2^2^\Q\ and so on. 
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Example 4.6. Consider the Gardner equation, which was shown to have two Hamiltonian structures with 

& = D x , S = AuD x + 2u x + 3s{uu x + u 2 D x ) - D 3 X . 

Hence, the operator connecting our hierarchy of approximate Hamiltonian symmetries is 

m = S ■ = 4u + Zeu 2 + (2 + Zeu^D- 1 - D 2 X . 

Therefore, our results on approximate bi-Hamiltonian systems will provide ready-made proofs of the existence 
of infinitely many approximate conservation laws and approximate symmetries for the Gardner equation. 
Note that the Frechet derivative for the right-hand side of Gardner's equation is 

T) K = 6(1 + 2eu)u x + 6(u + eu 2 )D x - D%, 

and 

ffi t = (4 + Qeu)u t + (2u xt + 3eu t u x + Seuu^D^ 1 
= 12uu x (2 + 5su) — (4 + 6eu)u xxx 

+ {Quu xx {2 + beu) + I2u x {l + 5eu) - u xxxx (2 + 3eu) - 3eu x u xxx ^J D x x . 

Theorem 4.7. Let Qq = eu x . For each k ^ 0, the differential polynomial Qk = & k Qo is a total x- 
derivative, Qk = D x Rk, and hence we can recursively define Qk+i = £%Qk- Each Qk is the characteristic of 
an approximate symmetry of the Gardner equation. 

Proof. To prove this theorem, we apply the similar method applied in theorem 5.31 of [5]. 

We proceed by induction on fc, so suppose that Qk = & k Qo for some Rk e si '. From the form of the 
approximate recursion operator, 

Qk+i = e(AuQ k + 2u x D~ 1 Q k - D 2 x Q k ) =eD x (2uD x 1 Q k + 2 J D~ 1 ( M g fc ) - D x Q k 

If we can prove that for some differential polynomial Sk e s/, uQk = D x Sk, we will indeed have proved that 
Qk+i = D x Rk+x, where Rk+i is the above expression in brackets. Consequently, the induction step will be 
completed. 

To prove this fact, note that the formal adjoint of the approximate recursion operator e5% is 

= e(4w - 2D' 1 ■ u x - D 2 X ) =D x x eS£D x . 

We apply this in order to integrate the expression uQk, by parts, so 

Qk = vM k [£u x ] = u x ■ {e@*) k [u\ + D x A k 

for some differential function Ak e si ' . On the other hand, using a further integration by parts, for some 
Bk e s/ the following identity holds: 

u x ■ {e@*) k [u\ = u x ■ D- 1 sM[u x ] = u x ■ D x x Q k = -uQ k + D x B k 

Substituting into the previous identity, we conclude 

uQ k = D x S k , where S k = -(A k + B k ), 

which proves our claim. □ 

Definition 4.8. A pair of skew-adjoint q x q matrix of differential operators 31 and S is said to form an 
approximately Hamiltonian pair if every linear combination aS>+b£', a, b e R, is an approximate Hamiltonian 
operator. A system of perturbed evolution equations is an approximate bi-Hamiltonian system if it can be 
written in the form 

c)ii 

— = KAu, e] m mjft * £SJ% (4.5) 
dt 

where @, § form an approximately Hamiltonian pair, and and Mi are appropriate Hamiltonian func- 
tionals. 
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Lemma 4.9. If 3, $ are skew-adjoint operators, then they form an approximately Hamiltonian pair if and 
only if @, $ and 3 + S are all approximate Hamiltonian operators. 

Corollary 4.10. Let Si and /§ be Hamiltonian differential operators. Then ^, $ form an approximately 
Hamiltonian pair if and only if 

prv w (e^) + prv, fl (0*) = o(e p ), (4.6) 

where 

Q CJ = I J{(9 a @9} dx, @g = i J"{6> a S 9} dx 
are the functional bi- vectors representing the respective Poisson brackets. 

Example 4.11. Consider the approximate Hamiltonian operators S>, *§ associated with the Gardner equa- 
tion. We have 

""»• - s °' <*»>. 4 - s d > (t *<«*) <4 

a, J J a,J \/3 = l / J 

in the case of the second approximate Hamiltonian operator for the Gardner equation, we have 

Wv se {u) = ^! Wv ge (u 2 ) = 2u£9, 

prv (?e (e@) = prv^ {9 a 9 X } dx = o(e p ) 
trivially, by the properties of the wedge product, it is deduced that: 

prv @e (6^) = prv^g J{(2u+ yu 2 )# a 6 x + ^9 X a 9} dx 



J 1 



{(2 + 3£M)^ a a X } = o(e p ) 
Thus ^ and S form an approximately Hamiltonian pair. 

Definition 4.12. A differential operator *2i : sd r — » srf s is approximately degenerate if there is a nonzero 
differential operator $ : .a?" 5 — > «af such that $ ■ $ = o(e p ) 

Now, according to [5], we are in a situation to state the main theorem on approximate bi-Hamiltonian 
systems. 

Theorem 4.13. Let 

u t = Ki[u,e] * Mifi <^<5J^ 

be an approximate bi-Hamiltonian system of perturbed evolution equations. Assume that the operator Qi of 
the approximately Hamiltonian pair is approximate nondegenerate. Let 3% = <§ ■ be the corresponding 
approximate recursion operator, and let Kq sb Qi8Mq. Assume that for each n = 1,2,... we can recursively 
define 

K n sb MK n ^ n>l, 

meaning that for each n, K n ^\ lies in the image of Si. Then there exists a sequence of functionals 
, 3tf\ , j • ■ • such that 

(i) for each n > 1 the perturbed evolution equation 

u t = K n [u,e\ * &5^f n * ALjC_i (4.7) 
is an approximate bi-Hamiltonian system; 

(ii) the corresponding approximate evolutionary vector fields v n = vx n all mutually commute: 

[v„, v r „] = o(e p ), n, m > 0; 
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(iii) the approximate Hamiltonian functionals ,¥6^ are all in involution with respect to either Poisson bracket: 

{Jt? n M 9 = o(e p ) = {Jf n ,J4? n } s , n,m:>0, (4.8) 

and hence provide an infinite collection of approximate conservation laws for each of the approximate 
bi-Hamiltonian systems (14. 5ft . 

We have seen that given an approximate bi-Hamiltonian system, the operator 3% = $ ■ , when applied 
successively to the initial equation Kq = !2)SJi?o, produces an infinite sequence of approximate generalized 
symmetries of the original system (subject to the technical assumptions contained in Theorem 14.131) . It is 
still not clear that Si is a true approximate recursion operator for the system, in the sense that whenever vq 
is an approximate generalized symmetry, so is v^q. So far, we only know it for approximate symmetries with 
Q = K n for some n. In order to establish this more general result, we need a formula for the infinitesimal 
change of the approximate Hamiltonian operator itself under a Hamiltonian flow. 

Lemma 4.14. Let ut = K *< SdSJV be an approximate Hamiltonian system of perturbed evolution equations 
with corresponding vector field v# = v# . Then 

prv^(fT) ^T> K ■@ + 3'-T)* K 

Theorem 4.15. Let ut = K «s 3>8,^{ kk S'SMq be an approximate bi-Hamiltonian system of perturbed 
evolution equations. Then the operators 3%i = e l S ■ < I s$ p, are approximate recursion operators for 
the system. 

Judging from 8£f = o(e p ), when I # 0, this type of approximate recursion operators have less significance 
than 8%q. 

Example 4.16. The approximate recursion operators of the Gardner equation are 

Bg Q = g . gr- 1 =4u + 2u x D- x + ZeiuUxD- 1 + u 2 ) - D 2 X1 8% x = e(4u + 2u x D- Y - D 2 X ) 

and we can apply 8$q to the right-hand side of the Gardner equation to obtain the approximate symmetries. 
The first step in this recursion is the flow 

Ut $53%{ * @8J%2 u xxxxx — 10 u u xxx — 20 u x u xx + 30 u 2 u x + e \55u 3 u x — 39uu x u xx — 9u 2 u xxx — 12m^ , 

which is not approximately total derivative, so we can not re-apply the approximate recursion operator to 
get a meaningful approximate generalized symmetry. 
But if we set 

Ki[u,s] = Q 5 = eKi[u,e] = e(6uu x - u xxx ), Jt = = eJ#b, Jt{ = 3*5 = eJ#{, 

then we can apply 3£q successively to K\ in order to obtain the approximate symmetries. The first phase 
become 

Ut * SdJti * ^(5^2 8%qK\ w e(u xxxxx — 10 uu xxx — 20 u x u xx + 30 u 2 u x ) 

in which 

■M = = I Uu 2 xx -5u 2 u xx + 5u 4 )dx 

is another approximate conservation law. 
Now, for Ki = 81qKi we have 

u t * S8Jt 2 « @5Jt 3 as & K 2 « e(-u xxxxxxx + l^uu xxxxx + A2u x u xxxx ) 

-1-70 s(u xx u xxx u u xxx ~\~ 2ii u x Auu x u xx u x ~j 

where 

f u 2 

Jt° 3 = 7e \ (-^ + uu 2 xx + 5 M 2 u 2 + u 5j dx 
is a further approximate conservation law. 
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5 Concluding Remarks 



Sometimes, differential equations appearing in mathematical modelings are written with terms involving a 
small parameter which is known as the perturbed term. Taking into account the instability of the Lie point 
symmetries with respect to perturbation of coefficients of differential equations, the approximate (perturbed) 
symmetries for such equations are obtained. Different methods for computing the approximate symmetries 
of a system of differential equations are available in the literature [IJ [2] [5] . 

The approximate symmetry method proposed by Fushchich and Shtelen [5] is based on a perturbation of 
dependent variables. This method has so many advantages such as producing more approximate group- 
invariant solutions, consistence with the perturbation theory, solving singular perturbation problems [101 II 1J 
and close relationship with approximate homotopy symmetry method [12] . But despite above mentioned 
benefits, this procedure converts a perturbed evolution equation to an equivalent perturbed evolutionary 
system. In his case, obtaining the corresponding Hamiltonian formulation will be hard. Due to increase of the 
dimensions of Hamiltonian operators f^, computation of the approximate recursion operator 8% = § ■ 
is difficult. 

Since prolongation and Frechet derivative of vector fields are linear, both of the approximate symmetry 
methods can be extended on the Hamiltonian structures. But due to the significance of vector fields in 
Hamiltonian and bi-Hamiltonian systems, the approximate symmetry method proposed by Baikov, Gazizov 
and Ibragimov [2 [2] seems to be more consistent. 
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